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ABSTHACT 


Until  Recently,  few  alternatives  existed  for  the  use  of 
hydraulic  and  pneumatic  actuators  in  primary  flight  control 
applications.  With  the  advent  of  the  sanarium-o  bait  perma¬ 
nent  magnet  brushless  dc  motor,  consideration  oust  now  be 
given  to  the  utilization  of  an  electromechanical  actuator  in 
missiles  which  require  significant  maneuvering  capability 
and  hence,  greater  torques.  Tris  tuesis  investigates  the 
theory  and  techniques  of  pulse  widtn  modulated  speed  control 
of  brushless  ic  motors.  After  describing  basic  pulse  width 
modulation  concepts,  two  constant  velocity  control 
schemes  are  presented:  current  feedback  and  a  limit  cycle 
scheme.  By  calculating  the  motor  form  factor  (a  figure  of 
merit  for  power  losses  in  the  switching  transistors  wnicn 
comprise  the  ? W M  network),  the  relative  worth  of  each  sen erne 
is  then  evaluated.  An  in  depth  study  is  conducted  of  the 
limit  cycle  approach,  with  an  emphasis  on  the  power  loss 
reductions  obtained  through  the  reduction  of  :  he  velocity 
limit  settings. 
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I.  INTRODUCTION 


Untix  recently,  few  alternatives  existed  tor  the  use  o: 
hydraulic  or  pneumatic  actuators  in  primary  flignt  control 
applications.  However#  advances  malt  in  the  field  of  rare 
earth,  permanent  magnet  materials  and  in  nigh  power  semicon¬ 
ductor  transistor  technology  has  made  possitlc  tne  ust  of 
electromechanical  actuators  as  a  practical  alternative. 
Utilizing  rare  earth  (specifically,  samarium  cobalt) 
magnetic  material  within  a  Lrusnless  dc  motor  design 
provides  a  prime.  rover  in  flight  control  applications 
offering  super ior  performance  characteristics  over  electroh- 
ydrauiic  systems. 

Elimination  of  the  brush  type  commutation  scheme  witnir. 
the  dc  motor  provides  numerous  advantages:  1)  higher  rated, 
motor  speed  along  with,  a  reduction  in  weight  and  volume  ror 
a  given  horsepower,  2)  the  anility  to  use  r ernanent  magnet 
rotors  instead  of  a  rotating  armature  winding,  combined  with 
the  elimination  of  the  brush  assembly  translates  into 
design,  implementation  and  maintenance  simplifications,  3) 
since  there  are  no  brushes,  no  arcing  will  occur,  hence 
allowing  motor  operation  in  hazardous  environments,  and  4) 
improved  thermal  characte  ristics,  as  losses  (ohmic  and 
core),  which  arise  primarily  within  tne  stationary  portion 
cf  the  machine,  are  easily  dissipated  through  the  stator 
housing . 

Use  of  the  brushless  dc  machines  does  not  come  without 
cerain  disadvantages,  chief  among  them  being  tne  cost  and 
the  uncertain  availability  of  the  samarium  cohalt  material 
foe  rotor  construction.  As  an  aside,  a  recent  study  was 
conducted  in  which  the  performance  characteristics  of  both  a 
ferrite  an  1  a  samarium  cobalt  type  dc  motor  w  *_-r  e 
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fol lowing  chapter.  The  research  was  conducted  utilizing  a 
computer  model  of  a  3-phase,  4  pole  brushless  ac  motor, 
which  us  developed  Ly  Thomas  in  a  related  study  [Ref.  2]. 

ft  is  to  be  noted  that  this  research  p  ro ject  represents 
only  ere  part  of  an  ongoing  effort  to  accurately  simulate 
t..e  use  cx  a  orushless  dc  motor  as  an  electromagnetic  actu¬ 
ator  for  use  in  advance  1  missile  control  systems.  While 
elect  rcnagr.eti  c  actuators  have  previously  seen  use  in 
missile  ;ro]t:ts  such  as  H ASH  and  Condor,  these  actuators 
have  teen  too  large  and  have  had  too  slow  a  response  for 
r.i  j  u  t  o  r  c  u  e  applications  as  are  x  o  uxi  u  i  a  t  n  e  A  .11 A  A  '■! 
(Advanced  Medium  Range  Air  to  Air  dissiie)  project  [Ref.  3]. 
Tre  overall  aim  of  the  entire  research  program  is  to  help 
exploit  the  technological  opportunity  which  exists  within 
lire  earth  permanent  magnet  1c  macnines  in  the  role  or  elec¬ 


tion.  echa.iicai  actuators  for  use 
dap  1 i cx t io  ns. 


m 


I  J  1  U  t)  ( 


missile  control 


II.  PULSEWIDTH  MODULATION 


Pulsewidth  modulated  switching  aaplif iers  order  cor.sii- 
erable  advantages  ir.  the  control  of  dc  motors.  Btrore 
describing  the  techniques  used  in  simulating  the  pulsewidth 
modulated  control  of  the  computer  model  mentioned  previ¬ 
ously,  it  will  be  important  first  to  outline  the  basic  prin¬ 
ciples  and  characteristics  of  pulsewidth  modulation. 

A.  PULSEWIDTH  MODULATION  PRINCIPLES 

The  pulsewidth  modulation  scheme  utilizes  transistors  in 
the  switching  mode,  whereby  the  transistors  are  switcaed 
into  and  out  of  saturation.  This  switching  action  results 
in  the  minimization  cf  power  losses  in  the  transistors,  vita 
a  savings  realized  in  reduced  heat  sinking  requirements  and 
in  the  usage  of  less  expensive  power  transistors .  Since  the 
power  transistors  are  switched  on  and  off  at  a  frequency 
beyond  the  system  bandwidth,  the  motor  will  filter  the  high 
frequency  components  of  t ne  modulated  signal  and  respond 
only  to  the  signal's  low  frequency  components. 

Closing  a  feedback  loop  around  the  pulsewidth  modulated 
amplifier  results  in  the  amplifier  behaving  as  either  a 
current  or  a  voltage  source.  The  feedback  loop  allows  the 
motor  to  be  connected  with  the  amplifier  through  additional 
series  inductance,  resulting  in  a  smoothing  of  any  current 
ripple.  The  feedback  loop  also  allows  for  easy  current 
limiting,  merely  by  limiting  the  output  of  the  feedback 
summing  amplifier.  Finally,  a  feedback  xoop  provides  output 
short  circuit  protection,  as  the  output  current  is  deter¬ 
mined  by  the  input  voltage  without  regards  to  tne  output 
impedance. 
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E.  MODULATION  TECHNIQUES 

Thtrc  are  two  basic  methods  ror  ODtair.ii.-j  a  pulse widtn 
mod Jiattd  sijnai. 

1  •  I ne  Either  let;. od 

The  first  technique  requires  that  tne  input  signal 
(7.0),  be  added  to  a  high  frequency  sawtootn  signal  (also 
known  as  a  dither  signal).  After  summing,  tne  resultant 
signal  ( Y 1 )  ,  is  fed  into  a  relay  element.  The  relay  element 
converts  the  summed  signal  into  a  two  level  output  (Y2)  , 
wnich  then  switches  from  +V  to  -V  whenever  Y1  experiences  a 
zero  crossing,  as  shown  in  Figure  2.1. 

The  duty  cycle  (a)  of  the  output  signal  is  related 
to  the  input  signal  and  magnitude  (E)  of  the  sawtooth  signal 
by: 

a  =  ( E  +  X)/2E  (eqn  2.1) 

Utilizing  this  technigue  has  as  an  advantage  the  fact  that 
one  may  control  the  frequency  of  tne  surplied  sawtooth 
signal  without  any  changes  to  tne  motor  controller 
circuitry. 

2.  The  Limit  Cycle  dethod 

The  second  technique  for  producing  a  puisewidta 
modulated  signal  is  by  closing  a  feedback  loop  around  a  two 
level  switch.  The  feedback  signal  causes  the  system  to 
exhibit  high  frequency  limit  cycle  behavior.  Taking  a 
velocity  control  system  as  an  example,  if  the  velocity  error 
signal  (reference  signal  minus  the  actual  velocity)  were 
negative,  indicating  that  the  velocity  was  too  high,  the 
output  signal  from  the  switch  would  tnen  serve  to  drive  tr.e 
velocity  lower.  As  the  velocity  of  the  system  dropped  below 


Figure  2.1  Creation  of  a  PH  M  Signal 
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ti.t.  rerorence  revel,  tne  error  signal  woard  1 1.  e  r.  become 
positive,  causing  the  ou.ut  of  the  switch  to  change  states 
ana  consequent ly  force  tae  veiocit/  to  once  again  increase. 
In  this  iianner  the  switch  output  wouii  thus  exhibit  oscilla¬ 
tions,  describing  a  rulsetrain  whose  frequency  as  determine  1 
by  the  voltage  levels  at  which  tae  switch  operates  and  by 
tine  dynamics  in  the  feedback  path.  The  frequency  of  tne 
velocity  waveform  would  necessarily  be  the  same  as  that  of 
the  output  of  the  switch,  and  under  steady  state  conditions 
would  become  constant  anc  periodic.  The  puisewidth  of  the 
signal  at  the  output  of  the  switch  is  dependent  a p on 
specific  system  dynamics.  For  tae  case  of  a  d.c.  motor 
operating  under  load,  the  puisewidth  is  directly  related  to 
the  load  on  the  motor.  The  specific  details  of  a  limit 
cycle  velocity  control  scheme  is  examined  in  Chapter  3. 


C.  ANALYSIS  OF  THE  PUISEWIDTH  MODULATED  SIGNAL 

A  diagram  of  an  ideal  dc  motor  is  shown  in  figure  *.2. 
The  "freewheeling"  diode  { FWD)  serves  to  bypass  the  motor 
during  the  pulse  off  reriod,  allowing  the  armature  current 
to  circulate.  Figure  2.3  shows  tne  typical  steady  state 
current  and  voltage  reiationsnip  in  a  puisewidtr  control 
scenario. 

Eecause  the  supply  voltage  is  being  switched  at  frequen¬ 
cies  typically  on  the  order  of  mag  iltude  of  5  KHz,  it  is 
important  to  study  the  power  losses  within  the  motor  with 
tne  power  being  pulsed  on  and  off.  A  first  approximation  to 
the  evaluation  of  power  losses  due  to  heating  within  the 
armature  resistance  is  seen  in  equation  2.2. 


f:  d2.  ..  ) 


( e  q  n  2.2) 


Figure  2.2  Basic  DC  Motor  Schematic 

Vie  may  now  define  the  current  form  factor  (k)  as  the  ratio 
of  the  BMS  current  to  the  average  current: 

*  =  Irri:./Iave  (e^n  2*3) 

Developed  motor  torgue  is  directly  proportional  to  motor 
current  (eguation  2.4),  allowing  motor  losses  (under  PWM 
conditions)  to  he  described  as  in  eguation  2.5. 

T  ,■  =  K;  *I  ive  (egn  2.4) 

P1  =  B*k2*I2aVe  (egn  2.5) 

Substituting  eguation  2.4  into  eguation  2.5  demonstrates 
that  motor  losses  are  dependent  on  tne  current  form  factor 
arid  the  output  motor  torgue,  as  seen  below: 


P 


(E/K,  2)  *k  2*T  2 


(egn  2.6) 


Figure  2.3  PHM  Steady  State  Behavior 


The  motor  form  factor  has  a  large  influence  on  motor 
heating  and  hence,  power  losses.  Since  performance  in  speed 
control  systems  is  ofttn  limited  by  power  dissipation 

constraints,  it  is  important  to  determine  the  form  factor 

for  a  given  PWM  scheme.  The  relationship  between  motor 
armature  losses  and  form  factor  is  shown  in  Figure  2.4 

[Bef.  4]. 

In  order  to  determine  the  form  factor  for  a  given 
system,  one  must  first  be  able  to  determine  a  system's 

average  and  RKS  currents.  The  differential  equations 
describing  the  motor  action  for  the  basic  dc  motor  system 
(as  seen  in  Figure  2.2)  are  as  follows: 


Pulse  on: 


Ldi/d t  =  V 


HI 


K,  *w 


(eqn  2.7) 


Figure  2.4  Additional  Armature  Losses  vs.  Form  Factor 


where  tn  is  the  system  electrical  time  constant. 


I  =  (V  -  K,  *w)/R  (egr.  2.11) 

and 


I„  =  *w/F  (e  gr.  2.  12) 


The  average  and  E n S  currents  are  then: 


I 

a  v  e 


(aln  -  bl.a)  /£ 


(egr.  2.  13) 


T2 


rms 


=  (1/P)  (  (aln  2  ♦  blm2)  -  t,  (I  ,)  (IJ) 


J 


(egr.  ^.14) 


where 


a  =  t  *ln  ( (I  -  I  )/  (I  -  I  )  ) 


(egr.  2.15) 


b  =  t  *ln  ((I  ♦  I  )/{I  +1)) 

ri  o  :n  oo  m 


(e'4».  2.  16) 


I  =  (I,  -I.) 


(egn  2.  17) 


I.  =  (I  ♦  I  ) 

•  r:i 


(egr.  2.18) 


and  1/P  is  the  switching  freguency. 
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As  u  means  of  demcnstra  tin  j  toe  viahixity  of  ?*;'  control 
of  dc  motors,  simulations  wore  conducted  using  ixXea  duty 
cycle  tover  pulses  to  determine  the  relations nir  between 
form  factor  and  motor  load  torque.  A  computer  program  was 
written  to  analyze  the  output  motor  current  waveshape  for 
the  average  and  ras  currents  utilizing  relationships 
detailed  in  this  section.  The  program  used  to  compute  these 
currents  may  be  found  in  Appendix  C.  Figures  2.5  and  2.c 
snow  that  the  form  factor  rapidly  approaches  unity  as  the 
load  on  the  motor  is  increased,  indicating  that  the  rotor  is 
experiencing  only  slight  additional  losses  (ir.  terms  of 
percentages)  due  to  the  power  pulsing  effect  as  compared 
with  a  constant  voltage  supply  arrangement.  It  is  also 

pointed  out  that  the  motor  runs  more  efficiently  (iower  form 
factor)  at  higher  frequencies  for  a  giver,  load  torgue.  Inis 
is  due  to  the  fact  that  as  the  motor  is  pulsed  core 
frequently,  the  motor  speed  will  not  drop  off  as  fast  and 
hence,  the  energy  required  to  move  the  mass  of  the  rotor 
hack  to  its  steady  state  speed  will  not  he  as  great. 
However,  switching  losses  in  the  transistors  will  usually 
limit  the  switching  frequency  to  less  than  ten  KHz. 
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Figure  2.6  For*  Factor  vs.  Load  Torgue  (5  KHz) 


III.  PCILSEKIDTH  MODULATED  SPEED  CONTROL 


T  ho  abiiit  y  of  a  dc  motor  to  maintain  a  jive,  a  speed  wren 
a  load  toroue  is  applied  is  generally  referred  to  as  speed 
regulation.  Although  a  do  motor  by  itself  is  an  open  loop 
system,  the  prescence  of  the  Lack  electromoti ve  force  (temf) 
signal  serves  to  close  a  natural,  "built-in",  feedback  loop, 
as  shown  in  Figure  3.1  However,  because  the  dc  motor  is 
intrinsically  an  open  loop  system  with  relatively  constant 
power  input,  as  the  load  torgue  increases,  the  speed  will 
decrease,  and  hence,  no  speed  regulation  may  be  achieved. 
In  order  to  maintain  a  constant  speed,  the  input  power  must 
increase  with  the  applied  load. 


Figure  3. 1  Block  Diagram  of  a  Basic  DC  Motor 

Figure  3.2  illustrates  the  basic  principle  behind 
puisewiith  modulated  amplifiers.  Rhat  is  important  to  note 


is  that  as  the  load  or*  the  system  increases,  the  duty  cycle 
cr  pulsewidta  of  the  input  signal  also  increases.  Cr.e 
thought  then,  is  to  attempt  motor  speed  control  ty  deter¬ 
mining  which  motor  parameters  are  changing  relative  to 
varying  loads  and  to  then  make  pulsewidth  a  function  of  one 
or  more  of  those  parameters.  3ecause  motor  current  varies 
linearly  with  load  torque,  some  form  of  current  feedback 
appears  to  he  the  logical  selection  for  implementation  of  a 
speed  control  scheme.  A  specific  current  feedback  control 
technique  was  investigated,  as  well  as  a  limit  cycle  control 
method,  the  details  of  which  will  now  be  presented. 


(a)  no  load  voltage 


(b/  heavy  loed  voltage 


tme 


!  *  me 


Figure  3.2  Pulse  Width  as  a  Function  of  Load  Torgue 
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A.  CURRENT  FEEDBACK 

E-rior  to  utveioj  ir.g  any  specific  control  schemes, 
sta lies  wert  cor.  lac  tel  or  lotcr  behavior  relative  to  a  f  x  :<e  1 
pulse  width  tower  pulsed  input,  it  a  ir  eguc-ncy  or  5  K  Hr. 
Motor  speed  was  round,  as  1:.  t re  rase  or  constant  sup  .  ly 
voltage,  to  vary  linearly  with  load  torque  ,  as  snow;,  it 
Figure  3.3.  As  ar.  aside,  all  studies  were  conducted  with 
the  load  tcroues  ranging  i  ror  zero  to  eighty  ounce -inch  es, 
as  this  ran^e  represented  the  linear  ran ge  or  operation  for 
the  motor  modelled  in  the  Thomas  study  [Ref.  2].  The  rlot 
of  average  current  vs.  load  tor  ana  for  the  fixed  puisewidtn 
sin  u  lat  x  cr.s  are  shown  in  Fie  are  3.4. 

Studying  the  curves  found  in  Figures  3.3  and  3.4  led  to 
the  conclusion  that  a  scheme  for  speed  control  of  tr.e  motor 
could  he  developed  with  the  *ulsewidth  being  made  directly 
proportional  to  tne  average  motor  current.  Tne  basic  forr. 
of  the  motor  pulsewidth  was  decided  to  be  as  follows: 

?  *  —  lC:  +  j.  _ . *  K  (egn  3.1) 

Tr.e  motor  parameters  found  in  equation  3.1  may  be  defined  as 
follows : 

?W  =  the  input  puisewidtn  (duty  cycle; 

I  ^  =  average  motor  current 

DCF  =  a  dc  term  which  establishes  no  load  speed 

K  =  the  current  feedback  constant 

It  was  decided  that  the  motor  would  run  at  a  minimum  of 
SOU  duty  cycle  pulses  to  minimize  power  losses  which  would 
occur  at  smaller  duty  cycles  in  light  loan  conditions.  Vith 
a  frequency  set  at  5  KHz,  this  necessarily  fixed  no  load 
speed  at  approximately  1375  rpm.  Since  lar^e  current  tran¬ 
sients  could  be  expected  when  the  motor  was  started  or  w,,ez 
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AVERAGE  SPEED  YS.  LOAD  TOEQl'E 


LFCKXD 

r»o  ~  rn.sF,  dctv  cyo 

7;V  I’l  LSF  I) l  TV  CY( 
NO  inv.M 


0  0  10.0  uo.o  :ui.o  10.0  r,  o.o  oo.o  to.o  no. 


LOAD  TORQl’K  (()'/— IN) 


Figure  3.3  notor  Speed  ?s.  Load  Toryue  (fixed  pulse  width) 


AVKRAOK  MOTOR  (TRRKNT  (AMI'S) 


MOTOR  (T/RRKXT  MS.  LOAD  TORQUE 


Figure  3.4  Current  vs.  Load  Torgue  (fixed  pulse  width) 
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tne  load  condition  was  changed,  a  limiter  was  written  into 

the  simulation  program  to  set  the  maximum  pulsewidth  at  just 
less  than  the  10  OS  duty  cycle  point.  It  was  not  set  at  1 00 
due  to  difficulties  enc  our.  ter  e  1  with  the  simulation 
language.  Other  additions  to  the  basic  computer  coded 
included  a  limiter  to  prevent  pulse  widths  with  less  tnar.  SO* 
duty  cycles  as  well  as  a  current  limiter  to  prevent  negative 
currents.  The  latter  was  added  to  simulate  the  effect  or 
the  addition  of  the  "freewheeling"  diode  described  ir.  the 
preceding  chapter. 

To  set  the  no  load  speed  at  1375,  it  was  necessary  to 
obtain  parameters  that  would  establish  PW  egual  to  .5  at 
zero  load  torgue.  The  value  of  K  was  determined  from 
Figures  3.3  and  3.4  by  noting  that  tne  75%  duty  cycle  condi¬ 
tion,  at  approximately  1350  rpm,  occurred  at  a  load  tcrgue 
of  40.0  oz-ins,  which  also  corresponded  to  an  average  motor 
current  of  2.50  amperes.  Setting  PW  to  .75  and  DCF  to  .43 
in  egn.  3.1  led  to  a  K  value  of  .107.  A  value  of  DCF  of  .43 
results  in  an  approximately  50%  pulsewidth  modulated  signal, 
as  no  load  motor  current  is  approximately  .002  amps.  It  was 
felt  that  higher  speeds  could  then  be  achieved  by  increasing 
DCF,  as  the  basic  relationship  between  pulsewidth  and  speed 
appeared  linear  for  any  given  load  torgue. 

Extensive  simulations  were  conducted  utilizing  the  feed¬ 
back  control  relationship  as  shown  in  equation  3.2.  Oraphic 
results  for  these  simulation  trials  are  shown  in  Figure  3.5. 

?  W  =  .48  +  ( I  *  .107)  (egn  3.2) 

It  is  cleariy  evident  from  Figure  3.5  that  the  control 
scheme  utilized  performed  u nsa tisf a ct or illy  for  its  task  of 
maintaining  constant  speed  throughout  tne  given  range  of 
load  torques.  While  certain  variations  might  have  been 
expected  in  the  output  speed,  tne  results  iemonstratei 


wuicii  were  clea 


non  iinearit  its  in  motor  performance 
unsuitarle  for  its  given  application.  One  of  the  major 
problems  encountered  with  this  control  scheme  stemmed  fro::, 
the  fact  that  the  averse  current  was  used  as  a  feedback 
parameter,  rather  than  the  actual  motor  current  (the  ripple 
present  in  the  motor  current  was  deemed  to  be  too  high  to  bo¬ 
used  in  a  velocity  control  scheme  eased  on  current  feed¬ 
back)  .  The  average  motor  current  proved  to  be  unsatisfac¬ 
tory  for  the  given  task  for  two  reasons:  1)  average  motor 
current  requires  time  to  approach  the  system's  actual 
average  current  value  due  to  the  cnanges  to  the  current 
incurred  by  transients  such  as  are  caused  by  changes  to  the 
system's  steady  state  behavior  and  thus  adds  significantly 
to  the  motor's  settling  time  to  variances  in  load  or 
commanded  speed,  and  2)  apparent  nonline-ant  res  which  would 
appear  if  pulse  duty  cycle  were  plotted  as  a  function  of 
load  tongue,  which  runs  contrary  to  the  initial  assumptions 
upon  which  the  control  scheme  was  devised. 

B.  VELOCITY  LIflIT  TECHNIQUE 

The  velocity  limit  control  scheme  has  teen  developed 
under  the  assumption  that  a  pure  velocity  command  has  been 
issued  by  the  motor  control  logic.  It  is  recognized  that 
other  systems  might  issue  torgue  commands  to  the  motor  in 
response  to  a  generated  missile  fin  position  error  signal 
and  the  current  motor  speed.  Simulation  of  the  complete 
electromechanical  actuator  is  not  the  intent  of  tnis  thesis, 
and  hence,  will  not  be  attempted  here.  The  schematic 
diagram  for  the  network  used  to  implement  the  limit  control 
scheme  is  shown  in  Figure  3.6.  This  network  will  control 
the  notor  speed  in  such  a  way  that  if  speed  is  below  its 
commanded  level  (plus  a  pre-defined  tolerance),  the  power  to 
the  motor  is  turned  on.  If  the  motor  speed  rises  above  t:.is 


Figure  3.6  Scheaatic  Diagram  Of  Pulsevidth  Modulator 

Operation  of  the  system  in  the  limit  mode  is  relatively 
straightforward.  Referring  to  Figure  3.7,  if  the  motor 
speed  is  below  {VCQd  ♦  VIOL),  transistor  QM  is  switched  on,- 
allowing  motor  currert  to  flew.  When  the  speed  reaches 
(VC 0,1  ♦  VTCL),  is  switched  off,  which  then  induces  a 
large  voltage  across  the  inductor  terminals,  owing  to  a 
rapid  rate  of  change  in  the  inductor  current.  This  induced 
voltage  turns  on  diode  DM,  whicn  provides  a  path  for  the 
decaying  motor  current.  When  motor  speed  decays  past  (VCOM 
-  VTCL),  CM  is  switched  LacK  on  again.  Transistor  Q E  and 
diode  DO  ace  utilized  when  the  motor  is  operated  in  the 
reg  o  nt  ia  ti  vt  braking  mode  of  operation.  In  this  mode,  the- 
dc  motor  is  used  as  a  generator  as  the  inertia  of  the  rotor 


Figure  3.7  Motor  Velocity  Waveshape 


is  converted  into  electrical  energy  to  charge  the  system’s 
dc  voltage  supply  while  controlling  the  braking  of  the 
rotating  masses.  Again,  how  this  mode  is  utilized  is  a 
function  of  the  design  of  the  motor's  electromechanical 
actuator,  and  will  not  be  further  commented  upon.  The 
control  of  transistors  Q B  and  QM  can  be  accomplished 
utilizing  voltage  comparators,  where  one  input  is  a  voltage 
proportional  to  the  commanded  speed  plus  the  allowed  speed 
tolerance,  while  the  other  input  is  a  voltage  proportional 
to  the  actual  motor  speed.  Waat  becomes  important  to 
realize  now  is  that  the  pulsewidth  and  the  pulse  freguency 
are  both  variable,  and  will  be-  dependent  on  certain  dynamics 
of  the  system. 
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A  trial  simulation  was  conducted  usm^  the  lxr.it  control 
scheme.  The  speed  tolerance  was  sot  at  five  rpm,  and  the 
commanded  speed  was  set  to  1400  rom.  Table  I  contains  the 


results  of  the  simulation. 


?his  data  is  also  represented 


graphically  in  Figure  3.3.  Loading  of  the  motor  was  accom¬ 
plished  using  a  terminated  ramp  signal;  the  terminal  value 
of  the  ramr  is  t ne  desired  motor  loading.  The  speed  accu¬ 
racy  of  the  motor  is  defined  as  the  difference  between  the 
maximum  and  minimum  motor  speed  divided  by  the  commanded 
speed . 

TABLE  I 


Motor  Characteristics 


Load  Torque  Speed  (rpm) 

(oz-in)  maximum  minimum 


0.0 

16.0 

32.0 

43.0 

64.0 

90.0 


1416.0 
1413. 2 
1417.  9 
1416  .  1 
1  4 1  4 .  3 
1409. 6 


1365.3 

1375.9 

1 309 . 7 
1366.2 
1362.  1 

1355.8 


'A  accuracy 


2.19 

3.00 

3.44 

3.57 

3.71 

3.79 


The  results  of  the  initial  simulations  indicates  that 


positive  control  of  the  commanded  motor  speed  may  be 
accomplished  utilizing  the  limit  cycle  method.  All  further 
studies  are  therefore  based  on  a  computer  model  whose  speed 
is  controlled  in  this  manner. 

Studies  of  system  cna racteristics  (pulsewidtn,  pulse 
frequency,  speed  regulation)  are  presented  in  detail  i  r.  the 
next  cnapter  in  order  to  tetter  define  the  operational 
enevelo^e  of  the  modelled  dc  motor  using  the  velocity  limit 
control  scheme.  Additionally,  tne  effects  of  the  addition 
of  series  inductance  arc  also  investigated  to  determine  its 
aufects  in  regards  to  ripple  reduction. 
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17.  PERFORMANCE  OPTIMIZATION 

The  simulations  described  at  tne  cr.  i  of  tne  previous 
cuapter  were  not  necessarily  indicative  of  the  o:  tiau;:. 
performance  potential  of  the  modelled  dc  motor.  It  appears 
evident  that  there  must  exist  means  by  which  the  current 
ripple  and  the  speed  regulation  may  be  respectively  reduce! 
and  improved.  One  way  to  achieve  improved  speed  accuracy 
(wnich  implies  that  the  ripple  present  in  the  motor  velocity 
is  reduced)  is  to  decrease  the  speed  tolerances  established 
for  motor  operation.  The  erfects  of  the  variation  of  tie 
speed  tclearr.ee  settings  will  be  examined  shortly.  However, 
to  reduce  the  current  ripple,  wnich  in  turn  translates  to 
reduced  power  losses  in  the  motor,  it  has  beer,  supjeste! 
tnat  cne  must  add  inductance  in  series  with  the  motor  arma¬ 
ture  [Ref.  5].  The  effect  of  the  addition  of  series  induc¬ 
tance  cr,  motor  performance  is  studied  in  the  f oil ewinj 
sec  t icn. 

A.  ADDITION  OF  SERIES  INDUCTANCE 

A  dimensionless  current  ripple  may  be  defined  as  the 
curer.t  ripple  multiplied  by  the  motor  armature  resistance 
and  divided  by  the  supply  voltaye.  The  current  rip  pit. 
itself  is  defined  as  the  difference  between  tne  current  at 
the  time  the  motor  is  pulsed  on  and  the  time  when  the  motor 
is  pulsed  off,  or  simply  as  the  difference  between  tne 
minimum  and  maximum  currents,  as  follows: 

I  =  i  (al)  -  i(0)  (e :  n  4.1) 
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g  =  1  -  exp (- t.j  (e an  4. a) 

Dimensionless  current  ripple  may  be  plotted  versus  tie 
puisewidth,  or  duty  cycle,  with  the  ratio  of  the  period  cf 
the  PWM  signal  to  the  electrical  time  constant,  tau,  to  :ors 
a  family  of  curves,  as  in  Figure  4.1. 

As  can  be  seen  from  Figure  4.1,  the  magnitude  of  the 
current  ripple  depends  to  a  great  extent  upon  the  ratio  of 
the  pulse  period  to  the  motor  electrical  time  constant.  V’e 
should  therefore  expect  that  that  reduction  of  the  motcr 
current  ripple  can  be  acco mplished  through  the  addition  of 
series  inductance,  which  reduces  tne  magintude  of  tau. 

The  effect  of  adding  series  inductance  to  the  motor  may 
he  seen  by  writing  the  differential  eguations  for  the 
simplified  circuit  diagram  of  a  dc  motor  as  shown  in  Figure 
4.2.  Applying  Kirchoff’s  Law  and  summing  voltages  around 
the  1  co  L.  results  in  a  cur  r  ent -vo  it  a  g  e  relationship  as  shown 
in  eg  uat  ion  4.5: 

V  =  ;  L  +  L 1 )  d  i/d  t  ♦  R  i  *  K .  •  w  ( e  j  n  -* .  5  J 

Taxing  the  Laplace  transform  of  equation  4.5  yields: 
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Figure  4 . 1 


Current  Hippie  vs.  Pulse  Duty  Cycle 


3  7 


Solving  t-jUatioi:  4.6  for  ao  tor  current  y  it  Ids : 

■.  r  ;■ ) ) .  c.  i  ;  -  .  :  ■  (.  ■ )  <eln  4-7> 

Fi  jt  equation  4.7  the  electrical  tine  constant  rs: 

t  .  =  ( L  +  LI)  /R  (e-;n  4  .  S) 


1 1-  e  t'*.  I  tCt:i 


t  he 


a  i  ii  t  io*.  or 


Strips 


1 1. u  ac- 


-O  test  ti.  e  e^rects  or 
tal.Ce,  tit  COu.t  itt  r  rOCC-i  WdP  J.OCificd  tj  Tell 
than  or  iGiactar.ce  e-ial  ir.  la^i.it  j.ie  to  the  mluctance 
present  ir.  the  tutor's  winJir.ys  (i  jr.  orin-j  the  cnar.ce  ir.  the 
total  resistance  for  tr.e  reasons  sentur.ea  aoeve)  . 
S;  ju  laticr.s  were  then  conducted  thro  a^hout  the  load  tor  .ue 


o : 


:ro  e 


ti'.'.ntv  ojnce-inci.es 


wi  t :. 


:no  tor 


ran  y  e 

velocity  set  at  at  1 403  rpa  (appr axina tel y  halt  of  the 
motor's  no  load  speed)  and  speed  tolerance  set  to  ±1  r  „.o . 
Table  II  summarizes  ti.e  results  of  these  simulation  trials. 
The  data  for  simulations  aa  de  without  t..e  additional  ir.iuc- 
tance  is  shown  in  Table  III  and  is  included  for  sake  of 
coij-arisor.  with  the  data  in  Table  II. 


TABLE  II 

Inductance  Effects  on  Motor  Operation 


Loa  u 

Ter  yue 

lave 

Ins 

Form  Factor 

Current 

(c  r  - 1 

ns) 

(aap-) 

(aap  s) 

Pipple  ( 

05.  3 

3.257 

3.4  54 

1 .  7b  5 

.1100 

16.  0 

3  .  995 

1.131 

1.  138 

.  0  9  74 

32.  0 

1  .  967 

2.  0  32 

1.  033 

.0^46 

43.  0 

3  .  0  23 

3.  0fc4 

1.  012 

.  0562 

64.  0 

4.039 

4.023 

1.005 

.  0737 

80.  0 

5  .  C  1 1 

5.  0  31 

1.  004 

.  1295 

TABLE  III 

Motor  Performance  (no  series  inductance) 


Lea 

(on 

1  Tor  out- 

"its)  “ 

Ta  vt 
( i  Q  L.  s) 

I  r.~is 
(a np  s) 

Fora  Factor 

Current 
Ripple  ( T) 

06. 

■3 

0.  179 

0.34  7 

1.933 

.32b 

In. 

0 

1.013 

1.10  1 

1.081 

.  255 

«i  • 

0 

1.363 

2.026 

1.022 

.  275 

*4  8  . 

J 

3.  0  1  3 

3.370 

1.003 

.  2  o  1 

t4. 

3 

4 .  0  2d 

4.040 

1.003 

.  2  2 1 

8  v . 

-j 

5.017 

5 .  0  2  6 

1.002 

.235 
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Figure  4.3  Inductance  Effects  on  Current  Ripple 
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Examination  of  the 


Tatles  II  an  1  III  lemon- 


data  in 

strates  clearly  the  ripple  reduction  obtained  through  the 
addition  of  the  series  inductance.  lotor  form  factor  was 
calcuiattvi  to  show  that  tower  losses  (at  low  load)  ir.  the 
transistors  are  also  reduced,  and  t.rus  represents  an  addi¬ 
tional  benefit  gained  from  the  added  inductance.  However, 
while  series  inductance  will  reduce  the  current  ripple 
witnout  affecting  steady  state  behavior,  it  will  have  an 
affect  on  the  transient  behavior  of  the  system,  as  will  no¬ 
te  shown. 

Assuming  now  that  the  system  is  operating  in  steady 
state,  a  step  input  command  (such  as  a  change  in  the  motor's 
commanded  velocity)  will  force  the  pud.sewidth  modulated 
signal  to  the  full  on  condition.  Ine  response  of  the  system 
will  now  be  limited  by  the  motor's  electrical  time  constant. 
Ir.  the  case  where  series  inductance  has  been  added  tc  the 
system,  the  response  of  the  motor  to  a  step  input  will  be 
slower  than  if  tne  additional  series  inductance  were  not 
pres  ent . 

To  test  this  effect,  simulation  trials  were  conducted 
where  the  motor  was  allowed  to  achieve  a  steady  state  speed 
of  1C00  rpm  and  then  was  subjected  a  stet  input  command  to 
increase  motor  speed  to  140  0  rpm.  Tor-ue  ioa  i  was  32  oz-ins 
and  the  speed  tolerance  setting  was  ±  1.  The  response  time, 
or  the  time  required  for  the  motor  to  settle  at  the  new 
commanded  speed,  was  measured  for  trials  in  the  motor's 
standard  configuration  and  for  the  case  where  series  induc¬ 
tance  was  added.  For  tne  first  case,  where  there  was  no 
additional  inductance,  the  response  time  was  measurea  at 
approximately  2.5  milliseconds.  When  series  inductance  was 
added  to  the  system,  the  response  time  slowed  to  approxi¬ 
mately  five  mi  11  is ec ends.  Thus  if  one  decid  es  to  add  series 
inductance  to  the  motor  to  reduce  tne  current  ripple 
effects,  that  decision  must  be  tempered  by  tne  fact  that  the 


transient  response  c£  the  system  will  ci.ange  due  to  the 
change  in  the  motor’s  electrical  time  constant. 

B.  F. EDUCTION  OF  VELOCITY  RIPPLE 

Optimization  of  motor  performance  will  require  taat  the 
ripple  content  of  the  motor  velocity  at  steady  state  he 
minimized.  To  reduce  the  ripple,  it  will  then  be  necessary 
to  reduce  the  speed  tolerance  settings  vmch  will  establish 
the  pulsewidth  modulated  input  signal  to  the  motor.  To 
determine  tne  effects,  if  any,  on  motor  performance  {other 
tnan  the  reduction  of  velocity  ripple),  it  was  necessary  to 
perform  a  number  of  simulations  with  varied  speed  tolerance 
set  tings . 

Three  simulation  trials  were  conducted,  with  the  speed 
tolerance  set  at  ±  5  rpm,  ±1  rpm  and  ±.  1  rpm.  Table  IV 
summarizes  the  results  of  the  simulations.  The  data  for 
Table  IV  is  plotted  in  Figures  4.4  and  4.5. 

TABLE  17 

Performance  Trials  for  Various  Speed  Tolerance  Settings 


Speed  Tolerance  (rpm) 

±5  ±1  ± .  1 


Load  Torgue 

rWle 

form 

rWle 

form 

r^fla 

f  orm 

(oz-in) 

factor 

factor 

factor 

05.  0 

2.  1  9 

2.  360 

0.39 

1. 765 

.214 

1.709 

16.  0 

3.00 

1.319 

1  .  16 

1.  138 

.242 

1.121 

32.  0 

3 . 44 

1.372 

1.22 

1.  033 

.235 

1  .026 

48.  0 

3.  57 

1. 283 

1.21 

1.0  12 

.228 

1.0  11 

64.0 

3.71 

1.095 

1.26 

1.005 

.228 

1 .0  05 

80.  0 

3.7  9 

1.054 

1.29 

1.003 

.250 

1.002 

It  is  apparent  that  reducing  tne  s^eed  tolerance  dees 
indeed  reduce  the  ripple  content  of  the  motor’s  speed. 
Additionally,  as  the  tolerance  is  reduced,  so  does  the  fore 
factor  of  the  motor(£or  specific  loading),  indicating  nigher 
motor  performance  efficiency.  Thus  as  speculated  before. 
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wit*,  lower  spet'J  tolerances,  tnere  exists  improved  motor 
perforaar.cc.  Of  course  this  should  nave  beer,  intuitively 
otvious  evt-r.  prior  to  coi.  ducting  tne  simulation  trials  due 
to  the  fact  that  the  greater  the  s^eed  tolerances,  the 
greater  the  inertia  obtained  by  tne  rotating  nass  of  the 
rotor  nefore  a  speed  limit  is  reached  ar.u  thus  greater 
amounts  of  energy  would  have  to  be  expended  in  order  to 
increase  the  motor’s  speed  back  to  the  current  commanded 
spee d . 

C.  CCMHENT  S 

’■'e  have  seen  where  the  simulated  performance  of  a 
modelled  brushless  dc  motor  has  been  improved  through  the 
addition  of  series  inductance  and  the  optimizing  of  the 
speed  tolerance  settings  of  the  chopper  control  logic. 
There  are  reasons  to  assume  that  the  performance  enhance¬ 
ments  noted  in  this  cnapter  may  not  necessarily  be  realized 
in  an  operational  cruise  missile  scenario.  The  final 
chapter  of  this  report  discusses  where  caution  need  be  taker, 
when  reviewing  this  work  and  prior  to  applying  these  results 
in  the  use  cf  brushless  dc  motors  in  aerospace  applications. 


Figure  4.5  Form  Factor  vs.  Load  Torgue  (VTOL  varied) 
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V.  RECOMMENDATIONS  for  further  studies 


As  stated  in  Chapter  3,  the  assumption  was  made  tnat  tno 
missiie's  ccntrol  logic  would  generate  velocity  commands  to 
the  motor  controller  unit.  This  nay  not  necessarily  re  true 
for  all  missile  applications.  Certain  controllers  nay 
generate  torque  commands  (or  equivalently,  current  commands, 
as  current  is  directly  proportional  to  motor  torque)  to  the 
motor  to  provide  a  constant  torque  to  affect  the  re  3 aired 
missile  maneuvers.  The  effect  that  this  would  nave  or.  the 
simple  ccntrol  scheme  studied  in  this  report  could  prove  to 
he  fairly  significant,  and  will  now  be  looked  at  in  closer 
detail. 

A.  BEYOND  SPEED  CONTROL 

One  major  change  that  would  oe  required  of  the  motor 
controller  if  torque  commands  are  to  be  issued  iron  the 
missile  control  logic  is  that  both  position  and  current 
loops  would  necessarily  have  to  be  closed  around  the  P'd'l 
amplifier.  Of  course,  the  implication  of  generating  tor3ue 
commands  is  that  the  operating  envelope  of  the  motor  by 
necessity  would  have  to  be  defined  to  include  the  plugging, 
braking  and  the  regenerative  braking  modes  of  operation. 
The  importance  inherent  within  the  inclusion  and  modelling 
of  these  nodes  is  test  justified  by  noting  that  the  regener¬ 
ative  braking  mode  serves  not  just  to  control  the  bracing 
torque  applied  to  the  rotor  but  also  to  recnarje  the 
missile’s  supply  batteries.  The  model  would  then  need  to 
have  logic  blocks  which  could  recognize  when  each  operating 
mode  would  be  appropriate  and  then  generate  tne  requisite 
commands  to  the  motor  control  logic.  The  idealized  step 


response  of  the  fin  actuator  system  to  comiaui  :hanct  in  zin 
position  is  shown  ir.  Figure  5.1  with  the  repaired  rotor 
operating  modes  shown  as  a  function  of  flap  position. 

Since  1)  a  position  loop  has  yet  to  be  closed  and  2)  no 
accounting  has  been  trade  for  the  modes  of  operation  ether 
than  motoring  in  the  forward  direction,  it  is  most  higr.Iy 
recommended  that  these  areas  be  investigated  to  more  accu¬ 
rately  model  the  motor  as  an  integral  part  of  an  overall 
missile  fin  actuation  system.  It  is  emphasized  here  that 
modelling  the  position  control  of  the  fin-motor  system  would 
represent  the  next  logical  step  ir.  accurately  simulating  the 
dynamics  of  the  electromechanical  actuator  sytem. 

B.  SIMULATION  DEFICIENCIES 

Studies  of  the  modelled  motor  nave  teen  concentrated  on 
the  response  of  the  system  to  stet  in^ut  commands  under 
constant  load  conditions.  Of  course,  under  normal  operating 
conditions,  the  system  most  undoubtanly  would  be  subject  to 
a  series  of  varied  load  conditions  as  the  missile  steers  it.: 
course  towards  the  target.  The  load  on  the  motor  would  then 
be  a  function  of  the  aerodynamics  to  which  the  missile  is 
subjected  during  flight,  i.e.,  missile  speed,  attitude, 
acceleration,  etc.  In  order  then  to  model  the  system  tetter 
under  tne  dynamics  or  flight,  the  effort  to  close  the  posi¬ 
tion  loop  should  be  followed  by  more  precise  load  studies, 
so  that  the  motor's  behavior  may  be  studied  within  a  context 
more  closely  related  to  its  predicted  operating  environment. 
Other  studies  that  would  prove  worthwhile  include  (but 
certainly  are  not  limited  to)  :  the  study  of  the  voltage 
switching  affects  on  the  controller's  power  transistors, 
designing  the  hardware  required  to  implement  the  motor 
controller  logic,  as  well  as  the  associated  software. 
Studying  of  the  effects  of  closing  a  pnase-lockel  servo  ioop 


F  f  n  Angle 

Regenerative  Braking  Mode 


Figure  5.  1  Theoretical  Fin  Actuator  Eesponse 
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aiSO  be 


considered  due  to  tee 


for  velocity  control  nijat 
system’s  precise  speed  regulation  capability. 

One  caution  worth  noting  is  that  the  model  used  as  a 
basis  for  all  simulation  studies  assumed  a  linear,  average 
flux  tack  erf  waveform.  The  back  emf  signal  is  in  actuality 
sinusoidal  in  nature,  and  would  therefore  result  in  the 
addition  of  fundamental  and  harmonic  frequency  components  to 
the  motor  parameters  which  were  studied,  suen  as  current  and 
velocity.  It  is  recommended  that  for  further  studies  a  mere 
auvanced  motor  model  which  simulates  sinusoidal  back  emf  be 
utilized. 

C.  SUMMARY  OF  RESULTS 

Pulsewidth  modulation  has  been  stown  to  be  a  viable 
method  of  accurately  and  reliably  controlling  tr.e  velocity 
of  a  brushless  dc  motor.  Form  factor  studies  indicated  that 
power  losses  in  the  switching  transistors  may  be  minimized, 
thus  allowing  for  smaller  power  transistors  and  reduced  heat 
sinking,  which  translates  into  reduced  costs  in  controller 
design  and  implementation.  The  addition  of  series  induc¬ 
tance  was  shown  to  a  have  a  very  definite  impact  cr.  motor 
current  ripple,  reducing  it  significantly  in  comparison  to 
simulations  conducted  without  the  series  inductance.  Jse  of 
the  limit  cycle  pulsewidth  modulation  scheme  was  shown  to  be 
a  superior  method  for  implementing  pulsewidth  modulation. 
Areas  where  further  research  efforts  may  continue  have  teen 
presented,  including  specific  recommendations  for  fcllow-on 
studies  to  this  thesis. 
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APPENDIX  a 
THE  MOTOR  MODEL 


The  motor  that  was  mole  lied  was  a  commercially  available 
trashless  dc  motor.  The  current  and  steed  carves  for  the 
motor  have  teen  shown  in  Figures  3.3  and  3.4.  Tre  motor  is 
a  three  phase,  four  pole  machine  with  the  commutation  being 
accomplisheu  electronically  utilizing  a  set  of  three  Hail 
effect  position  sensor  devices  and  a  set  of  six  switching 
transistors.  The  switcaing  takes  place  every  30  degrees  of 
mechanical  rotation. 

The  tack  EMF  signal  was  assumed  to  be  directly  propor¬ 
tional  to  the  motor  speed.  The  actual,  sinusoidal  nature  of 
the  waveshape  was  not  taken  into  account  in  the  model  used. 
What  follows  is  a  brief  description  of  certain  procedures 
that  were  added  to  the  basic  motor  simulation  program,  a 
description  of  program  variables  added  to  the  initial 
program  as  well  as  observations  made  concerning  the  execu¬ 
tion  of  the  simulation  program. 

A.  PROCEDURES  ADDED  TO  THE  BASIC  PROGRAM 

1.  2£22®3ure  IC  LIP 

Procedure  ICLI?  was  added  to  the  model  to  account 
for  the  addition  of  the  freewheeling  diode  FWD.  Currents 
were  then  necessarily  clipped  (negatively  going)  at  zero 
amperes,  thus  preventing  the  circulation  or  negative 
currents. 

2.  Procedure  VC1P 

This  procedure  takes  as  an  input  the  velocity  error 
signal  and  yields  the  motor’s  input  voltage  based  u^ou 


COu  a  n  d  U  S  .  t't'  0.  and  t  f.e  dDlXS  i.ti  S  pee  a  tolerance.  at  i 

with  in  this  ;r  octdurt  that  the  limit  cy  cl  e  Lt:.avior  or  th 
sys  tea  —  s  e  s  t  a  b  a.  rs  n  e  d . 

3  .  Proctiure  RESET 

As  the  output  of  the  integrator  oioc,<  which  yield 
rotor  position  counts  upward  continuously  from  zero, 
procedure  was  required  to  reset  the  rotor  position  to  z er 
decrees  when  360  degrees  of  rotation  of  tne  rotor  wa 
achieved.  RESET  keeps  track  of  tne  nurarer  of  tines  tnat  th 
rotor  turns  past  the  360  degree  point  and  uses  this  informa 
tion  to  update  the  variable  TiiRST,  which  cycles  ranges  fro 
zero  to  360  degrees.  The  importance  of  this  procedure  wil 
ultimately  he  realized  when  a  position  loop  is  closed  wit hi 
the  system. 

B.  MOTOR  PARAMETERS 

The  parameters  which  follow  are  those  added  to  the  tasi: 
program  to  achieve  the  required  speed  control  effects. 

1 .  VCMD 

VCMD  is  the  commanied  motor  velocity. 

2 •  1IQL 

VIOL  is  the  velocity  tolerance  limit  in  rpm. 

3 .  ICL  IP 

ICLIP  is  the  motor  current  which  has  oeen  ad  juste 
to  prevent  negative  current  flow. 

4 .  THR  ST 

THRST  is  the  rotor  position  in  degrees  which  range, 
from  zero  to  363  degrees. 


c 


NCTES  ON  PROGRAM  EXECUTION 


A  fixed  interval  integration  technique  was  used  in  lieu 
of  the  variable  step  Runge-Kutta  method  supplied  as  a 
default  integration  technique  in  C3MP-III.  The  trarc?oicai 
technique  was  used  as  it  demonstrated  itself  to  perform  as 
accurately  as  tne  varuable  step  methods  but  used  better  than 
507  less  computer  time.  The  integration  interval  was  cnosen 
as  .000001  seconds. 

Since  the  switching  frequency  of  the  motor  was  on  the 
order  of  2500  KHZ,  to  accurately  observe  motor  behavior 
during  the  pulse  on  and  off  periods  a  print  interva-  of  less 
than  50  microseconds  (typically  20  microseconds  was  chosen) 
was  required.  Because  of  this,  it  was  difficult  to  observe 
the  micr cscopic  detail  of  motor  operation  in  terms  of  the 
variances  present  in  the  current  and  speed  for  periods  of 
greater  tnan  one  second,  as  the  OSIP  program  is  limited  to 
approximately  5500  lines  of  output.  For  the  studies  male 
for  this  report  the  limitation  encountered  did  not  pose  a 
major  problem,  but  could  prevent  an  obstacle  to  further 
studies.  Of  course,  as  tr.a  studies  of  this  system  advances, 
the  requirements  for  such  detailed  assessment  of  motor  oper¬ 
ation  may  not  be  r resent,  and  larger  print  intervals  may  be 
used,  thus  allowing  studies  to  be  of  greater  duration. 


//AGKCSM  P2  jGd  ( 1 102,  0  1 1  oi  t'ASKINASArtAA*  »CLASS=C 

//♦MAIN  0RG=NPGVM  1.  llOIr 

//  EXEC  CSMFXV 

//X. COMPRINT  OC  ClMMY 

//X.SYSPRINT  uC  CuMMY 

//X.SYSIN  CO  * 

INJTI AL 

CONSTANT  nT  =  15.9,  oM  =  0.00015,  El  =  0.0,  c l  =  0.0, 

M  -  1.0,  JM  =  C.C01,  KB  =  0.112,  T LP  =  10a. 

PARAMETcR  LA  *  .0016,  RA  =  2.740,  ... 

TLP  =  05.,  VTCl  =  1.,  VLMC  =  1400. 

♦  A1  =  LA/RA  --  T PE  INVcfiSE  ELECTRICAL  TIME  CONSTANT 

♦  A2  =  J/b  —  IhE  INVEkSE  MECHANICAL  TIME  CONSTANT 

♦  VTuL  —  THE  V  EluC  IT  r  LIMIT  GETTING 

♦  VCMD  —  THE  CLmMANOEG  motor  velocity 

♦  TLP  --  lCAG  tcrgue  in  02 -ins 

♦  KB  --  EACK  EMf  CONSTANT 

♦  KT  --  TuRCoE  CONSTANT 

♦  N  —  GEAR  R  A  T  10 

♦  JM  --  PGTCR  MOMENT  jf  INERTIA 

♦  TlP  ~  LCAO  TCRuoE  w hcRE  PEAK  P0*ER  CcT  CCCJRS 


NO  SORT 

dLP  =  BL/(N**2J 
JLP  =  JL  /  ( N*  * 2  J 
J  =  JM  +  JLP 
B  =  EM  +  EL  P 
A 1  =  2.0  *  LA  /  RA 
A2  =  J  /  B 
THRST  =  0.0 
JF  AC  =  0  .C 

DYNAMIC 

VERR  =  VC  MO  -  mMRFM 
V I F=  vCLP  *  S1EPIC.OJ 
TL1  =  RAMPIO.CJ 
TL2  =  RAMP(.Cl) 

TL  =  T  L^  *  (Til  -  TL  2  J  /.  Cl 

VIB  =  0.0 

VI  N  =  VI  F  *  vie 

VlNl  =  VI N  -  VEMF 

VIN2  =  VIM  *  11. 0/RAJ 

IM  =  REALFL(0.CtAl,VIN2» 

TM  =  ICL  IP  *  K  7 

♦  TNI  =  TM  -  TL 

TN  2  =  TN  1  *  (  1 .0/BJ 

wM  =  REAL  PL  iQ.3»A2,T  N2J 

RMRP  M  *  *M  *  (  iO.  /3.  1415  9  ) 

wMRPMR  =  RMkPM/N 

VEMF  =  wM  ♦  KE 

ThETA  =  IMuHlC.C.^M) 

TnuEG  =  THETA  *  (  laO  .  0 /3  .  14  1  5  92c  5  J 


TnCON  =  7hR  S  7 
PwR  =  wM  *  7  M 


44  44  4  4  4  4  44  4  44  44  *  4  4  4  44  4  4  $***■*  *  $*  $  44  4 

*  PROCEDURE  ICLIP  SIMULATES  PkESlENCE  OF  4 

*  FREERFEE  L  INC  CICCE  dY  CLIfPING  T  *-  E  NEGATIVE  * 

*  ouING  CURRENT  A7  C  AMPS  * 

44  4  444  4  4  44  4  444  4  <  4  4  4  44  « 

PROCEDLRE  ICLIP  =  CLPCRIIM) 

IF  (IM  .  L  £  .  C.C)  ILLI  P  =  0.0 
IF  4 1  M  •  G  7  .  0.0)  ICLIP  =  IM 
ENDPRGCEDURE 


4444  44  4  4  44  444444444  4«=  $$****:*$$#  *¥*#«:;!;(=<:**:(;*■¥¥*$$=!'**** 

*  PROCEDURE  VCeF  ESTABLISHED  LIMIT  ^YCLE  dEFAVlCR. 

*  IF  THE  VELOCITY  IS  dE LOw  The  ESTABLlSFEO  SETPOINT, 

*  THE  INPuT  VOLTAGE  IS  PulSED  ON  (  VCLP  *  30)  ARC  IF 

*  THE  VELOCITY  IS  AEjV  E  VT  Cl  +  VCMC,  THE  INPuT 

*  VOLTAGE  IS  PLLSEO  OFF  (VCLP  =  0.0) 

44 4444 44444 44 44 4444 444444444444 444444 4 44444444 4 4444 44 


PRGCEDIRE  VCLP  =  WV(vERft»  VTCL) 

MVTCL  =  -1.  C  *  V7CIL 
IF  (VERR  .  G  7 .  VTCL)  VCLP  =  30.0 
IF  (VERR  .LT.  MVTGL)  VCLP  =  0.0 
ENDPRGCEOuRE 

4444444444 444 4  4  4  4 4 4  444  44 44  4 44 4  4 4444444444 4  4*444444 

*  TNI  SETS  THE  SIuN  OF  THE  LGAC  TORwUE  TG  ENSURE  * 

*  THAT  ThE  TGRCLE  OPPOSES  ThE  MGTOR  TQKCuE  RhE  N  * 

*  THE  MOTOR  OPERATES  IN  THE  FORWARD  DIRECTION  * 

*  AND  ACCS  TG  7 F E  MOTOR  TORCUE  *HEN  * 

*  THE  MOTOR  OPERATES  IN  THE  REVERSE  DIRECTION  * 
$**£****^***£*****4,**  4444444 444444444 4444 44444444 


PROCEOLRE  TM  =  f  FCet*C(  VIN.TM.TL) 

IF(  V  IN.L  T.O.C  )  GO  TO  10 
TNI  =  TM  -  7L 

GO  TO  1  £ 

10  TNI  =  TM  +  TL 
15  CONTINUE 
ENDPRGCEDURE 

♦*#***<**#  *=f>$  **$$•$**$*  4444 4  44 44444444 4444444444  44  4 

*  PwR  A  HAS  EE  EN  ADDED  TO  THE  PROGRAM  Tb  ACCOUNT  * 

*  FOR  ThE  NONLINEAR  ASPECT  OF  THE  MODELLEC  * 

*  MOTOR'S  OCT  PL  7  PORER  FUNCTION  * 

44 4444  4444 444444444 44 44444444 4  44444444 444444444  44  4 


PROCEDURE  PhR  A=  7  wIDCU  PRR, 7  L,  TLPi 
IP(TL.GE.TLP)  GC  TO  20 
P  w R A  =  PwR 
GO  TO  25 

20  TwlC  =  4.  50  *  (TL  -  TLPi 
PwREX  =  EXP(-TRIU) 

PwRA  =  PwR  *  PwREX 


####« 


2  5  CuM  IN  IE 
ENCPkGCEDURE 


*  ThIS  PRGCEDuRE  *AS  ADDED  Tu  RESET  THE  RGTCR'S  * 

*  POSITION  AFTER  IT  REACHES  3c0  DEGREES  * 

*  BACK  TG  0  DEGREES  * 

PkGCEDCRE  THR  S  T=R  E  S  £T  (  JFAC  »  THDEuJ 
TS  =  JFAC  *  3EC.0 
THRST  =  ThDEG-  TS 
IFITHRST  .LT.2EC.0J  GO  TO  AD 
JF  AC  =  JFAC  +  1  .0 
40  CONTINUE 
cNOPRCCEGuRE 


TERMINAL 

TITLE  BASIC  CC  MCTuR  SYSTEM 

TIMER  FINTIM  =  .015,  CuTCEL  =  .00004,  ... 

PRDEl  =  . OOOC A »  DECT  =  .  000001 

METHCD  TRAP 2 

*  Output  *mrpm,  theta,  ThRST 
PRINT  ICLIPtbUP.  WMRPM,  THRST 
LABEL  MGTCR  SPEEC  Du E  TO  STEP  INPuT 
PAGE  MERGE 

*  PAGE  XYPLOT 

END 

STOP 

ENDJ  GB 
/* 
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$E  NT  R  Y 

$J06  CuRR  »XR  EF 

C 

C  THIS  PROuRAM  IS  DESIGNED  TO  YIELD  Tnt  AV  ERAoE  CURRENT 
C  AND  THE  RMS  CURRENT  AS  1ft  EL  L  AS  K— F  a  U  TGk  FOR  GIVEN 
C  SIMULATION  Of  A  COMPUTER  MODEL  OF  A  P«M 
C  CONTROLLED  6  PC  ShLESS  DC  MCTGR 
C 

FEAL  IC,  ICC,  RES,  INDJC, RPMAVE, Pi, TIME, TIMES, FRtu 
REAL  EG.  IAv»I2kMS,A»d  ,ZN  ,KFAC  ,  I  \MS  ,  ICN,ICFF,IM»IN 
REAL  VIN  ,7C  ,P  ,T  L  »CC.  Pm»  ICEN 
PI  =  3.1415S 

RES  IS  TFE  RESISTANCE  CF  Th£  COILS 
INGUC  IS  THE  INDUCTANCE  EETaEEN  TFiE  TwO  TERMINALS 
VIN  IS  TFE  INPUT  VGL  TAoE  CF  THE  MOTuR 


RES  =  3.3AC 
INDu  C  =  .CO  13 
VIN  =  30. C 

FREC  =  FRECUENCY  OF  PwM  *AVE  SHAPES 

OC  REFERS  TC  CuIV  CYCLE  CR  *  CN  TIME*  OF  THE  WAVE 

TL  IS  TFE  L C AC  TCRUUE  FOR  TFE  GIVEN  RUN 

10  IS  TFE  PEAK  CURRENT  ACHIEVED 

IOC  IS  TFE  M  IN  IX  UM  CURRENT  SEEN  GY  THE  MOTOR 

TU  IS  THE  AMCLNT  OF  TIME  ThAT  TFE  CURRENT  WAVESHAPE 

IS  LESS  THAN  ZERO 

FREw  =  1163. 

CC  =  .654 
TL  =  S.6 
RPMAVE  -  14C0 
IC  =  A.A1 
ICC  =  3.42 
TO  =  .OOC 

£N  IS  THE  SYSTEM  ELECTRICAL  TIME  CONSTANT 

i N  *  INLLC/RES 
F W  =  (1.0/FRECJ  *  DC/100. 

TIME  =  PW 

TIME  2  -  1 . C/FREU  -  PR 
EO  IS  THE  BACK  EMF 

EC  =  1.631  *  RPMAVE  *.0558  *PI/3C.O 


IN  IS  THE  THEORETICAL  MAXIMUM  OF  THE  CURRENT 
IM  IS  TmE  THEORETICAL  MINIMUM  uF  THE  CURRENT 
A  IS  THE  PULSE  CN  TI ME 
B  IS  THE  PULSE  CFF  T IME 


IM  =  EC/FES 

IN  =  (VI  b  -  ECJ  /RES 

A  =  ZN  *  A  L  C  G  ( (  IN  -  ICu)/(IN  -  ILJi 
E  =  Z N  *  ALCGillQ  *  If'i/llJD  *  IMil 


F  -  A  ♦  E  ♦  TC 
I A  V  =  < (A*IN)  -  (E*IM JJ  /  P 
IDEK  =  ZMIIC  -  I  GO  )  *  ( I N  ♦  I  Mi  i  /  P 
I2RMS  =  ((A*IN*IN  +  B*IM*IMJ  -  I  DEN 
IRMS  =  SCF7(I2RMS) 

KF  AC  =  IFME/IAV 


RR  IT  E  (  6  ,1C  J 

10  FCRMAT(1X  .'STATISTICS  FCR  P«M  CONTROLLED  CC  MOTOR:'! 
FRITE  (6,iCiFPEW 

20  FCfiMAT(  IX, 'FREQUENCY  =  •  ,  F8 . 2  , '  H  Z  •  , ) 

RRITE  (  6  ,22  i  T  L 

22  FORMAT(  IX, 'LOAD  TGRQLE  =  ',F8.2i 
RR  IT  E  (b  ,50)  IN,  IM 

50  FORMAT  (IX, 'IN:  *,F10.3,'  IM:  »,F  10.3) 

RRITE  (  6,  ICC)  IAV.IRMS 

100  FORMAT  (1  X  j*  I AV  =  ',Fa.4,/f'  IRMS  =  *,Fb.4i 
RRITE  ( 6,2C0)  KF  AC 

200  FORMAT  (IX,  4K  FACTOR  =  ',F6. 51 
STOP 
EMC 

SENTRY 
$E  NTRY 


TIME 

ICLIP 

VCLP 

i^RP  P 

.3 

.0 

20.000 

.0 

4. CC300D-05 

. 3£6o  0 

30.000 

1.12  20 

8.0C0C0C-05 

.7*44  6 

20.000 

4.4261 

1.2C300C-04 

1 .  C6  7  7 

30.000 

5.86  56 

1.6CCC0C-04 

1.3484 

30.  000 

17.  224 

2.0 C000  l-0  4 

1 . 71o8 

20.000 

26. 765 

2.4CC00D-04 

2 . 02  3  0 

30.000 

38.  066 

2.8CG00C-04 

2.31 71 

30.000 

51.221 

3.2C000C-04 

2.3993 

20. OCO 

66.059 

3.6C000D-04 

2  •  6o  9  7 

30.000 

82.647 

4.0C300C-04 

2.  1285 

30.000 

100. 79 

4.4COOOC-04 

3.2759 

30.000 

120 .47 

4. 8C0C0C-04 

2 • cl  1 9 

30.000 

141.61 

5 . 2  COOO  C-0  4 

2.836  9 

30. OCO 

164.  13 

5.6CCOOD-04 

4.C510 

30.000 

187 .55 

6 . 0 COOO  C-0  4 

4 .2544 

30.000 

213.10 

6  .4C0OOC-O4 

4.447  1 

20. OCO 

239.40 

6. 8COOOC-04 

4.6295 

30.000 

266 . £4 

7.2C00OC-O4 

4.6018 

20.000 

295 . 25 

7.6COOOC-04 

4.9641 

30.000 

2  24 . £6 

8.  0  COCO  C-0  4 

5.1165 

30.000 

255 . 22 

8 .4  COOO C-0  4 

5.2594 

20. OCO 

2  8o . 68 

8.6COCOD-04 

5.392  9 

30.000 

418 .60 

9 . 2  COCO  C-0  4 

5.5172 

2C.000 

451.62 

9.6  0000  C-0  4 

5.662  5 

30.000 

485 .49 

1.0  C0CQD-03 

5.7391 

30.000 

519.64 

1 .04000C-03 

5.6370 

20. OCO 

554. 79 

1 . 0  EOOO  C-0  3 

5.526  6 

30.000 

590.21 

1. 120G0C-03 

6.C08  0 

30. OCO 

626 .24 

1.160 00  C-03 

6 .0814 

30. OCO 

662 . 63 

1.2  COCOG-Q3 

c. 147  1 

30.000 

699 . 74 

1.240 00 C-0 3 

6.2052 

30.000 

737. Cl 

1.2  €000  C-0  3 

6  .2560 

30.000 

774.61 

1. 32000C-03 

6 .2996 

30.000 

812  •  48 

1 .3  6000  C-0  3 

6.33o 2 

20. OCO 

650. 59 

1. 4COCOD-03 

6.3062 

30.000 

663 .89 

l . 4  4000  C-0  2 

6.2895 

30.000 

527 .25 

1  .  4  €  0  00  C  -0  3 

6 .406  6 

20. OCO 

565 .53 

1. 520C00-03 

6.4175 

30.000 

100  4  .o 

1 • 560CQC-03 

6.422  5 

20.000 

1043.3 

1. 6C0C0  C-03 

6.421  7 

30.000 

108  2  .0 

1.6 4CC0C-03 

6.4154 

30.000 

1120.6 

1 .6 6000 C-0 3 

6.403  7 

30. OCO 

1155.2 

1. 720C00-03 

6.3868 

30.000 

1197.7 

1.7  6000C-02 

6.2650 

20.000 

1236  .1 

1.8  COOO  C-0  3 

6.2384 

3C.0C0 

1274  .4 

1. 640000-03 

6.3072 

30.000 

1312.4 

1 . 8  £000  C-03 

6.2716 

20.000 

135 C  .2 

1. 92000 C-03 

6.231  7 

30.000 

1387  .8 

1. 960C0C-03 

5.5513 

.0 

1424.7 

2  .  OCOOO  C-03 

5.543  1 

.0 

1455.2 

2 • 0  40C00-03 

5. 143  8 

.0 

149 1.3 

2.0  6000  C-03 

4 . 753  4 

.0 

152 C  .9 

2. 12000  C-0  3 

4.2722 

.0 

1546  .2 

S3 


2.  U0C0t-03 

4 . COO  0 

.0 

1573.2 

2.2C00QC-03 

3.6369 

.0 

1596  .0 

2 . 2  40 COD-O 3 

3.2829 

.0 

lol  6  .5 

2 . 2  60C0  C-03 

2.9380 

.0 

1635.0 

2.320C0C-03 

2.602  2 

.0 

lc5  1 .3 

2.360C0  2-03 

2.2755 

.0 

lbo  5  .7 

2  .4C000C-03 

1.9379 

.0 

167  6.1 

2 . 44GCQD-03 

1.649  3 

.0 

1 6c  8  .5 

2.460 002-03 

1.349  7 

.0 

lc9  7  .2 

2.520CQC-03 

1 .059  1 

.0 

1709  .0 

2. 560000-03 

.  7  775  1 

.0 

1709 .1 

2.6CCC0C-03 

•  5C47  4 

.0 

1712.4 

2. 64000  C-0j 

.24080 

.0 

1719.2 

2.6  60C0C-03 

.0 

.0 

1719.3 

2 .72000  C-0  3 

.0 

.0 

1713.8 

2. 7  6CC0D-03 

.0 

.0 

1713.3 

2.6C0OOC-O3 

..c 

.0 

1712.7 

2. 84000 C-0 3 

.0 

.0 

1712  .2 

2.8  60C00-0  3 

.0 

.0 

171  l.o 

2.920COC-03 

.0 

.0 

1711.1 

2. 960COC-03 

.0 

.0 

171C.5 

3.0CCOOC-03 

.0 

.0 

1709  .9 

3 .04000  C-0 3 

.0 

.0 

1709.3 

3.0  6CC00-03 

.0 

.0 

1706.7 

3 •  1 2000  C-0  3 

.0 

.0 

170  6.1 

3. 16000 C-0 3 

.0 

.0 

1707  .3 

3 . 2  CO  000-03 

.0 

.0 

1706.9 

3.24OC0C-03 

.c 

.0 

1706.3 

3.260C0C-03 

.0 

.0 

170  5  .6 

J.32QC0C-03 

.0 

.0 

1705.0 

3  .3  6000  C-0  3 

.0 

.0 

1709.3 

3 • 4  CO  COD-03 

.0 

.0 

1703.7 

3.44000C-03 

.c 

.0 

1703.0 

3. 460  00  C-0  3 

.0 

.0 

170  2  .3 

3. 52COOD-03 

.0 

.0 

170  l.o 

3 . 56CC0  C-03 

.0 

.0 

1701.0 

3. 6CCG0C-03 

.0 

.0 

1700.3 

3.640COC-03 

.0 

.0 

1699.6 

3  .o  6000  C-0  3 

.0 

.0 

1696.8 

3.  7  20000-03 

.0 

.0 

1698.1 

3.7  60000-03 

.0 

.0 

1697.4 

3.8COOOC-03 

.0 

.0 

1696  .7 

3.  84CC0Q-03 

.0 

.0 

1695.9 

3.860C0C-03 

.0 

.0 

1695.2 

3. 92000C-03 

.0 

.0 

1694.4 

3.96GC0C-03 

.0 

.0 

1693.6 

4. 0  0000  C-0  3 

.0 

.0 

1692.9 

4.0  40  COD-03 

.0 

.0 

lo9  2.1 

4.0  6000  C-03 

.0 

.0 

169  1.3 

4. 12000C-03 

.0 

.0 

1690.5 

4.  1  6000  C -0  3 

.0 

.0 

lo8  9  .7 

4 . 2C0C0  C-03 

.0 

.0 

1688.9 

4.240C0C-03 

.0 

.0 

1686.1 

4.2  60C0C-02 

.0 

.0 

168  7.3 

4. 320  00  C-0  3 

.0 

.0 

1686  .4 

4. 36CC0C-03 

.0 

.0 

1685  .6 

4.4C0C0C-03 

.0 

.0 

1684.7 

4.44000  C-03 

.0 

.0 

I6tt3  .9 

4.460000-03 

.0 

.0 

168  3  .0 

4.520C0C-02 

.0 

.0 

1682.1 

4.5  6GC0C-03 

.0 

.0 

led  1 .3 

4.6CGCOC-02 

.0 

.0 

168  C  .4 

4  •  b  4  0  00  C-03 

.0 

.0 

1679  .5 

4.  6  60C00-03 

.0 

.0 

lo7  8  .6 

4.7  2300C-03 

.0 

.0 

1677.7 

4. 7 6000 C-0 3 

.0 

.0 

1676.7 

4.8  C0000-03 

.0 

.0 

lo7  5  .8 

'I 


mr--i w  ~  w— ■» 
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• 

4. 8  4000  C-Oj 

.0 

.0 

1674  .9 

4.8  6GC0C-03 

.0 

.0 

1673  .9 

4.9  2000  C-0  3 

.0 

.0 

1673.0 

4. 960  COD-03 

.0 

.0 

1672.0 

c 

5.0CGG0C-03 

.0 

.0 

1671  .1 

5  . 04000  C-0  3 

.0 

.0 

1670 .1 

5.06CCQQ-03 

.0 

.0 

16o  9.1 

5. 120COC-03 

.,0 

.0 

1666.1 

5.  1 6CC0C-03 

.0 

.0 

166  7  .1 

5 . 2  C  C  CO  C-0  3 

.0 

.0 

lcc  6.1 

3  .2  4000  C-0  3 

.0 

.0 

1665.1 

5.2  4  C  COD-O  3 

.0 

.0 

loo  4 .1 

A 

5.3  2Q00D-03 

.0 

.0 

16o3.1 

5 .3  6000  C-0 3 

.0 

.0 

166  2.1 

. 

5.4  CC00D-03 

.0 

.0 

loo  1  .0 

5.44CC0C-03 

..0 

.0 

lco 0.0 

r 

5.4  6CC0C-03 

.0 

.0 

1656  .9 

5. 520C0C-03 

.0 

.0 

1657.  y 

5  .56  000  C  -0  3 

.0 

.0 

165  6  .6 

V 

5. 6  C  C  COQ-03 

.0 

.0 

lc5  5  .7 

0( 

5.O4CC0C-03 

.0 

.0 

1654.6 

6.O60C0C-03 

.0 

.0 

1653.5 

5. 7  2C00D-0  3 

.0 

.0 

1652  .4 

3. 760  CO C-0 3 

..0 

.0 

165  1.3 

5.  8CCC0C-03 

.0 

.0 

1650.2 

5.6  4CG0C-03 

.0 

.0 

1649.1 

5. 660  00  C-0  3 

.0 

.0 

1647.9 

5. 92000C-03 

.0 

.0 

1646  .o 

• 

5.96CC0C-03 

.0 

.0 

1645.7 

6.  CC0C0C-03 

.0 

.0 

1644.5 

b.  C4GC0C-03 

.0 

.0 

1643.3 

6 . 06000  C-0  3 

.0 

.0 

1642.2 

6. 120COO-03 

.0 

.0 

1641.0 

b.l  6 OCO  C~ 0  3 

.0 

.0 

1639.8 

6. 200  00  C-0  3 

.0 

.0 

1638  .6 

O.240C0D-03 

.0 

.0 

163  7.4 

Q 

6.26000C-03 

.0 

.0 

1630.2 

6.  320G0C-03 

.0 

.0 

16^>  5  .0 

b. 36000C— 03 

.0 

.0 

lo33.o 

6.4  0000  C-0  3 

.0 

.0 

1632.5 

6. 440 CO D -03 

.0 

.0 

163 1  .3 

O.46000C-03 

.0 

.0 

1630  .1 

6. 5  2000  C-0  3 

.0 

.0 

1626.8 

6.560000-03 

.0 

.0 

162  7.5 

H 

6 . 6  COCO  C-0  3 

.a 

.0 

1626.3 

6.64000C-03 

.0 

.0 

1625  .0 

6  •  6  £CC0o-03 

.0 

.0 

lc2  3.7 

6.7  2000  C-0 3 

.0 

.0 

1622.4 

6. 7  6 CCOD-03 

.0 

.0 

1621.1 

6.8C000C-03 

.0 

.0 

1619.8 

o  .84000C-03 

.0 

.0 

1616  .5 

6.860000-03 

.0 

.0 

1617.2 

• 

6 . 9  2000  C-03 

.0 

.0 

1615.6 

6.9  6000  C -03 

.0 

.0 

1614  .5 

7 . 0  CCG0C-Q3 

.0 

.0 

1613.2 

7 .04000  C-0  3 

.0 

.0 

161 1  .8 

7.060000-03 

.0 

.0 

1610.4 

7.12000  C-0  3 

.0 

.0 

1609 .1 

7. 16000  C-0  3 

.0 

.0 

160  7  .7 

► 

7 . 2  C0C00-O3 

.0 

.0 

1606  .3 

• 

7 . 2  4000  C-03 

.0 

.0 

lc04 .9 

7.2  6000C-03 

.0 

.0 

1603  .5 

7.32000C-03 

.0 

.0 

1602  .1 

7.3c000C-0  3 

.0 

.0 

160  C  .7 

hi 

7.400000-03 

.0 

.0 

1599  .3 

7 . 4  4000  C-0  3 

.0 

.0 

1597.9 

0 

7.4  6000  C-0  3 

.0 

.0 

1596  .4 

CO 

> 

► 

4  \  1 

--  --  •  - 

•  -  i  ^ _ -■  -  ^  . -  V.  .'j 

7 . 5  2000C-0  3 
7. 56000C-03 

7.6  CGOOD-O  3 
7. 64CC0  C-03 
7.68CC0C-03 

7. 7  * COOu-Od 
7 .76000C-03 
7. SCO 00 C -03 
7 • 6  4000  C-0  3 

7.  88CC0C-03 
7 . 52CC0C-03 
7.96000  C-03 
8 . 0  COCOD-03 
8.040C0C-03 
8 .06000  C-03 
8. 1 20  C0D-03 
8. 163C0C-03 
d. 2  COOO  C-03 
8.2  4CC0C-02 

8. 2  8000  C-0  3 
8 . 320000-03 

8.3  60C0C-03 

8  .4  COOO  C-0  3 
8. 440000-03 
8.480COC-03 

8.5  2000  C-03 

8.5  6C00C-03 
8.60000  C-03 

8.6  4CC0D-03 
8 . o  8000  C-03 
d.  7  2000  C-03 

8.7  60  COO-O  3 
8.80000C-03 

8.  8  4C00  C-03 
8.88000  C-0  3 
8. 9  2000  C-0  3 
8.96CC0D-03 

9  .00000  C-03 
9. 04000 C-0 3 
9.0  8000D-03 

9 . 1  2000  C-03 
9. 1 6000  C-03 

9.2  CGG0C-Q3 
9.2 4000  C-03 

9.2  80000-03 

9.3  2000  C-03 

9. 3  6300  C-0  3 
9.4C003Q-03 

9 . 4  4000  C-03 

9.4  8000C-03 
9. 520COC-03 

9.5  6000  C-0  3 
9. 6CCC00-03 

9 . 6  4000  C-03 
9. 680  00  C-0  3 

9.7  20C0D-03 
9. 760 CO C-0 3 
9. 8C000C-03 

9.8  4GC0C— 03 

9.8  8000  C-03 
9. S20C0D-03 

9.9  6000C-03 
i  .  000 00  C-0  2 
1. 004000-02 
1.0C8C0C-O2 
1.01200  C-02 
1.0  16C0C-02 


1595.0 
59  -  .5 

1592 .1 
159  C  .o 

1509.1 
1587  .o 

1586.1 
158  4  .o 

1583.1 
158 1  .6 

1580.1 
1578.0 
157  7.0 

57  5  .5 
574.0 
57  2.4 
57C.8 
56  5  .3 
567.7 
566 .1 
564  .5 
5o  2  .9 
56  1  .3 
155S.O 
1558.0 
1556  .4 

1554.7 

1553.1 

155 1.4 
1545.0 

1548.1 

1546.4 

1544.7 
1543.0 

1541.3 

1535.6 
153  7.9 

1536.2 

1534.4 

1532.7 
153  C  .9 

1525.2 
152  7.4 

1525.6 

1523.9 

1522.1 

1520.3 

1518.5 

1516.7 
'514.0 

512  .0 

1511.2 

1505.3 
150  7.5 
150  5.6 
1503  .8 

1501.9 
1500. 0 
1498.1 
1496  .2 
1494  .6 

1492.4 
1490  .5 
148  8  .6 

1486.7 

1484.7 
1482  .8 


1 • 0  2CC0D-02 

.0 

.0 

l4o  0  •'y 

1.02400C-02 

.0 

.0 

1479.0 

1. 026000-02 

.0 

.0 

147  7  .0 

1.C3200C-02 

.0 

.0 

1*75  .1 

1.03600  C-02 

.0 

.0 

1*73.2 

1 . 04CC0C-02 

.0 

.0 

147  1 .  j 

1.044C0C-02 

.0 

.0 

1469  .* 

1  .04800  0-0  2 

.0 

.0 

146  7  .4 

1. C  £2000-02 

.0 

.0 

l4o  5  .5 

1 .056000-02 

.0 

.0 

1*6  3  .6 

1.06000  0-0  2 

.0 

.0 

14o  1.7 

1. 064000-02 

.0 

.0 

1459  .7 

1  .068000-02 

.0 

.0 

1457.8 

1.072000-02 

.0 

.0 

1455.9 

1.076  CO  0-0 2 

.0 

.0 

1454  .0 

1  .08000  0-02 

.0 

.0 

1*52  .0 

1.0  84C0D-02 

.0 

.0 

1450.1 

1. 088000-02 

.0 

.0 

1448  .2 

1. 092000-02 

.0 

.0 

1*46  .3 

1. 096000-02 

.0 

.0 

14*  4  .4 

1 . 1  COCO  C-02 

.0 

.0 

1442.4 

1.104000-0^ 
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